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Abstract 

Factorization of an n x n unitary matrix as a product of n diagonal 
matrices containing only phases interlaced with n — 1 orthogonal matrices 
each one generated by a real vector as well as an explicit form for the Weyl 
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1 Introduction 



Matrix factorization is a live subject of linear algebra. It seems that no general 
theory is yet available although many results appear almost every day. However 
our goal will not be so ambitious to present a general theory of matrix factor- 
izations but to tackle the problem of factorization of unitary matrices. Unitary 
matrices are a first hand tool in solving many problems in mathematical and the- 
oretical physics and the diversity of the problems necessitates to keep improving 
it. In fact the matrix factorization is closely related to the parametrization of 
unitary matrices and the classical result by Murnagham on parametrization of 
the n- dimensional unitary group U{n) is the following: an arbitrary nxn unitary 
matrix is the product of a diagonal matrix containing n phases and of n{n — l)/2 
matrices whose main building block has the form 

jj _ ( cos 6 —sin6e~'^'^ \ 

\sine e'^ cos 9 J ^ ^ 

The parameters entering the parametrization are n{n — l)/2 angles 6i and n{n + 
l)/2 phases ipi. For more details see [Q. 

A selection of a specific set of angles and/or phases has no theoretical sig- 
nificance because all the choices are mathematically equivalent; however a clever 
choice may shed some light on important qualitative issues. This is the point of 
view of Harari and Leurer ^ who recommended a standard choice of the Cabbibo 
angles and Kobayashi-Maskawa phases for an arbitrary number of quark genera- 
tions and accordingly they propose a new parametrization (factorization) which, 
up to a columns permutation, is nothing else than Murnagham parametrization 
with the change 6 — >■ i[/2 — 9 and 99 + vr. 

However there are some problems that require a more elaborate factorization. 
To our knowledge one of the first such problems is that raised by Reck et al. 
who describe an experimental realization of any discrete unitary operator. Such 
devices will find practical applications in quantum cryptography and in quantum 
teleportation. Starting from Murnagham parametrization they show that any 
nxn unitary matrix An can be written as a product 

where Bn E U (n) is at its turn a product of n — 1 unitary matrices containing 
each one a block of the form (1) and C„_i is a f/(n — 1) matrix. Consequently the 
experimental realization of a n x n unitary operator is reduced to the realization of 
two unitary operators out of which one has a lesser dimension. The experimental 
realization of a U{3) matrix, sketched in Fig. 2 of their paper, suggests that it 
would be preferable that phases entering the parametrization should be factored 
out, the device becoming simpler and the phase shifters, their terminology for 
phases, being placed at the input and output ports respectively. 
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A mixing of the Murnagham factorization and that of Reck et al. is that 
proposed by Rowe et al. [§] in their study on the representations of Weyl group 
and of Wigner functions for SU{3). The last parametrization is also used by 
Nemoto in his attempt to develop generalized coherent states for SU (n) systems 

i- 

Another kind of factorization is that suggested by Chaturvedi and Mukunda 
in their paper aiming at obtaining a more "suitable" parametrization of the 
Kobayashi-Maskawa matrix. Although the proposed forms for n = 3,4 are aw- 
fully complicated by comparison with other parametrizations existing in litera- 
ture, and for this reason it cannot be extended easily to cases n > 5, the paper 
contains a novel idea namely that that an SU (n) matrix can be parametrized by 
a sequence of n — 1 complex vectors of dimensions 2,3, ... ,n. Fortunately there 
is an alternative construction as it may be inferred from the construction of an 
SU{3) matrix as a product of two matrices each of them generated by a three- 
and respectively two-dimensional complex vector in a much more simple form 
than that presented in ref. |0. 

One aim of this paper is to elaborate this alternative construction in order to 
obtain a parametrization of x n unitary matrices as a product of n diagonal 
matrices containing the phases and n — 1 orthogonal matrices, each of them 
generated by a n- dimensional real vector. As a byproduct we obtain the Weyl 
form IP of a unitary matrix W = w*dw where w is a unitary matrix, w* its 
adjoint and d a diagonal matrix containing phases. The Weyl factorization was 
the key ingredient in finding the "radial" part of the Laplace-Beltrami operator 
on U (n) and SU{n) ||, |10| and this explicit form could help in finding completely 
the Laplace-Beltrami operator on unitary groups. 

The paper is organized as follows: In Sect. 2 we derive a factorization of 
n X n unitary matrices as a product of n diagonal matrices interlaced with n — 1 
orthogonal matrices generated by real vectors of dimension 2, 3 . . . , — 1. The 
explicit form of orthogonal matrices is found in Sect. 3 and the paper ends with 
Concluding remarks. 



2 Factorization of unitary matrices 

The unitary group U{n) is the group of automorphisms of the Hilbert space 
(C", (■,■)) where (■, ■) is the Hermitian scalar product {x,y) = Yl^jZ^'^Vi- If 
An G U{n) by A* we will denote the adjoint matrix and then A^An = In, where 
In is the nx n unit matrix. It follows that det An = e"^, where if is a phase, and 
dirriR U (n) = n^. 

First of all we want to introduce some notations that will be useful in the 
following. The product of two unitary matrices being again a unitary matrix 
it follows that the multiplication of a row or a column by an arbitrary phase 
does not affect the unitarity property. Indeed the multiplication of the j^^ row 



3 



by e*'''^ is equivalent to the left multiplication by a diagonal matrix whose all 
diagonal entries but the j*^ one are equal to unity and ajj = e^'^K The first 
building blocks appearing in factorization of unitary matrices are diagonal ma- 
trices written in the form dn = (e**^^, . . . , e**^") with ipj G [0, 2 vr), j = 1, . . . ,n 
arbitrary phases, and all off-diagonal entries zero. We introduce also the nota- 
tion d^''' = (1„_A:, e^^^, . . . , e*'^'-'), k < n, where l„_fc means that the first {n — k) 
diagonal entries are equal to unity, i.e. it can be obtained from dn by making the 
first n — k phases zero . Multiplying at left by dn an arbitrary unitary matrix the 
first row will be multiplied by e*'''^, the second by e**^^, etc. and the last one by 
Multiplying at right with d^'^ the first n — k columns remain unmodified 



and the other ones are multiplied by e 



#1 



giVfe respectively. A consequence 



of this property is the following: the phases of the elements of an arbitrary row 
and/or column can be taken zero or vr and a convenient choice is to take the 
elements of first column non-negative numbers less than unity and those of the 
first row real numbers. This follows from the equivalence between the permuta- 
tion of the i^^ and j*^ rows (columns) with the left (right) multiplication by the 
unitary matrix Pij whose all diagonal entries but an and equal to unity, 

an = ajj = 0, a^j = a^j = 1, i j and all the other entries vanish. In conclusion 
an arbitrary An G U{n) can be written as a product of two matrices, the first one 
diagonal, in the form 

An = dn An (2) 

where An is a matrix with the first column entries non-negative numbers. 

Other building blocks that will appear in factorization of An are the rotations 
which operate in the i,i + 1 plane of the form 



J, 



( h-1 



\ 







COS Oi 
sin Oi 



—sin Oi 
cos Oi 












,n 



(3) 



The above formula contains the block (1) with phase zero unlike other parametriza- 
tions in our parametrization the phases will appear only in diagonal ma- 

trices. 

Let V be the vector v = (1, 0, . . . , 0)* G 5'2„,-i G C" where t denotes transpose 
and S2n-i is the unit sphere of the Hilbert space C"' whose real dimension is 2?t,— 1. 
By applying An & U (n) to the vector v we find 

/ an \ 



AnV 



where a G S2n-i because An is unitary. The vector a is completely determined 
by the first column of the matrix An- Conversely, given an arbitrary vector of 
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the unit sphere w G S2n-i this point determines a unique first row of a unitary- 
matrix which maps w to the vector v. Therefore U{n) acts transitively on S2n~i- 
The subgroup of U{n) which leaves v invariant is U{n — 1) on the last n — 1 
dimensions such that 

S2n-i = coset space U{n)/U{n — 1) 

A direct consequence of the last relation is that we expect that any element 
of U{n) should be uniquely specified by a pair of a vector h G S2n~i and of an 
arbitrary element of U{n — 1). Thus we are looking for a factorization of an 
arbitrary element An & U (n) in the form 

An = Bn-(l / ^ ] (4) 



n-l 



where Bn E U (n) is a unitary matrix whose first column is uniquely defined by a 
vector b G S2n-i, but otherwise still arbitrary and An-i is an arbitrary element of 
U{n — 1). For the SU{3) group such a factorization was obtained recently @, 0]. 
Iterating the previous equation we arrive at the conclusion that an element of 
U{n) can be written as a product of n unitary matrices 



where 



An = Bn ■ -Bn-l • • • -B" ^ (5) 

,k I Ik 



Bn-k 



Bn- 



k 



Bk, k = 1, . . . ,n—l ,are kxk unitary matrices whose first column is generated by 
vectors bk G S2k-i', for example is the diagonal matrix (1, . . . , 1, e*'''"("+i)). 

The still arbitrary columns of Bk will be chosen in such a way that we should 
obtain a simple form for the matrices -B^"'', and we require that Bk should 
be completely specified by the parameters entering the vector bk and nothing 
else. In the following we show that such a parametrization does exist and then 
An G Un in Eq.(5) will be written as a product of n x n unitary matrices each 
one parametrized by 2/c — 1, k = 1, . . . ,n, real parameters such that the number 
of independent parameters entering An will be 1 + 3 + ■ ■ ■ + 2n — 1 = as it 
should be. 

In other words our problem is to complete an n x n matrix whose its first 
column is given by a vector 6„ G S2n-i to a unitary matrix and we have to do it 
without introducing supplementary parameters. For n = 3 this was found by us 
in |jTl| in an other context and here we give the construction for arbitrary n. 

If we take into account the property (2) the problem simplifies a little since 
then 
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where the first column of Bn has non-negative entries. Denoting this column by 
bi we will use the parametrization 



bi = {cos 9i, cos 62 sin 6'i, . . . , sin 6'i . . . sin 6n-if (6) 

where 6'i G [0, 7r/2], i = 1, . . . , n— 1; we call 9i angles. Thus Bn will be parametrized 
by n phases and n — 1 angles. According to the above factorization Bn is nothing 
else than the orthogonal matrix generated by the vector bi. Thus with no loss of 
generality Bn = dn On with On G 0{n). In this way the factorization of An will 
be 

-An = dn On C^n-l • • • ^^2 ^2 ^1 

where C^.^ has the same structure as -B^.^, i.e 

Consequently the factorization of unitary matrices reduces to the parametrization 
of orthogonal matrices generated by an arbitrary vector of the real n-dimensional 
sphere and in the next section we show how to do it. 



3 Parametrization of orthogonal matrices 

An operator T applying the Hilbert space Ti in the Hilbert space Ti.' is a contrac- 
tion if for any v E H, ||Tt;||7^/ < ||f||>{, i.e. ||T|| < 1, [|12|. For any contraction 
we have T*T < I-^/ and TT* < ly^ and the defect operators 

Dt= [I-H-T^Tf^, Dt* = {Iw - T T*Y'^ 

are Hermitian operators in Ti and Ti' respectively. They have the property 

TDt = Dt* T, T* Dt* = Dt T* (7) 

In the following we are interested in a contraction of a special form, namely 
that generated by a n-dimensional real vector b G R", i.e. T = (61, ... , where 
bi are the coordinates of b; its norm is ||T|| = {b,b) and T will be a contraction 
iff (6, b) < 1, i.e. if 6 is a point inside the unit ball of R". If (6, b) = 1, that is the 
case we are interested in, T is an isometry, i.e. 

T*T = 1, and Dt = 

and in this case Dt* is an orthogonal projection. A direct calculation shows that 
det {\In — Dt*) = A(A — such that the eigenvalue A = has unit multiplicity 
and the eigenvalue A = 1 is degenerated. From the second relation (7) we have 

Dt* T = DT*h = T Dt = bDT = 
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i.e. b is the eigenvector of Dt* corresponding to A = eigenvalue. 

The orthogonal matrix On which brings the operator Dt* to a diagonal form 







is the orthogonal matrix we are looking for because it is generated by an arbitrary 
n-dimensional real vector of unit norm. The multiplicity of A = 1 eigenvalue being 
n — 1 the form of the matrix On is not uniquely defined. Indeed ii vi, . . . ,Vn are 

the eigenvectors of Dt* 



Dt*vi — DT*h — 0, DT*Vk — Vk, k — 2, . . . ,n 
orthogonal eigenvectors will be also the vectors 



(8) 



vi and Gk = Rvk, k 



2,3, 



n 



where R G 0{n) is an arbitrary rotation acting only on the last n—1 eigenvectors. 
Thus there is a continuum of solutions for the orthogonal basis that diagonalizes 
Dt* ■ In this situation we have to make a choice between the possible bases. Our 
criteria was that the resulting orthogonal matrix On should have as many as 
possible vanishing entries. We found such a matrix that have (n — l)(n — 2)/2 
null entries and the result is expressed by the following lemma: 

Lemma l:The eigenvectors of the eigenvalue problem Eq.(8) which are the 
columns of the orthogonal matrix On £ SO{n) generated by the vector parametrized 
by Eq.(6) are given by 



( cos 9i 
sin 9i cos 02 



V2 = 



/ —sinOi 
cos 9i cos 62 
cos 9i sin 62 cos 6*3 



Vk+2 



\ sin 01 . . . sin 9n-i ) 



( Ofc 

-sin 9k+i 
cos 9k+i cos 9k+2 



\ 



\ cos 9i sin 92 ■ ■ ■ sin 9n-i / 



(9) 



k^ l,...,n-2 



\ cos 9k+i sin 9k+2 ---sin 9n-i J 

where 0^ means that all the first k entries are zero and k — 1,2, . . .n — 2. Alter- 
natively 

d 



d9k 



vi{9. 



^fe_i = 7r/2) A; = l,...,n-1 
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The full 0{n) group is obtained by multiplying On given by Eqs.(9) with the 
diagonal matrix which has all the entries hut the last one 1 and dn,n — 

Proof. Elementary checking shows that (f vj) = 6ij, i, j = 1, . . . ,n, and thus 
Vk are hnearly independent. Because the multiphcity of the null eigenvalue is 
unity it follows that Vk, k — 2, . . . ,n, are orthogonal eigenvectors corresponding 
to A = 1 eigenvalue. 

Lemma 2 The orthogonal matrices On ( 0^_^) at their turn can be factored 
into a product of n — 1 (n-k-1) matrices of the form Ji,i+i; e.g. we have 

On — Jn-l,n <^n-2,n-l ■ ■ ■ Jl,2 (10) 

where Ji^i+i are n x n rotations introduced by Eq.(6). 

Remark. If the angles that parametrize On are ^i, . . . , 6'„_i, then the angles 

that parametrize are denoted e.g. by . . . , 92n-3, etc. and the last angle 

entering O^'^ will be 0n{n-i)/2- 

Putting together all the preceding information one obtains the following result 
Theorem: Any element An E U (n) can be factored into an ordered product 

of 2n — I matrices of the following form 

An = dn On Cl ^^-1 ■ ■ ■ dr^O^'d^' (10) 

where d'^_^ are diagonal matrices and 0^_^ orthogonal matrices whose columns 
are given by formulae like (9) generated by real {n — k)- dimensional unit vec- 
tors. Using factorization (10) the above formula can be written as a product of n 
diagonal matrices and of n{n — l)/2 rotations Jk,k+i- 

The condition Y^'^lH'^^^^'^ Vi = 0, imposed on (pi the arbitrary phases entering 
the parametrization of An, gives the factorization of SU(n) matrices. 
If Wn = On dl_, Ol_, . . . dr'Or'dr ' = a dl_, Wn-i then 

Wn^W^dnWn (11) 

is one (of the many possible) Weyl representation of unitary matrices. 

If all the phases entering An are zero = 0, i = 1, . . . , n{n + l)/2, one gets 
the factorization of rotations Rn G SO{n) 

Rn = OnO\_,...Or' (12) 

and the full group 0{n) is obtained by multiplying (12) with a diagonal matrix 
d = (1, . . . , 1, —1) that has one entry equal to —1. 

Remark. The above factorization is not unique and we propose it as the stan- 
dard (and simplest) representation. Equivalent factorizations (parametrizations) 
can be obtained by inserting matrices like Pij as factors in the formulae (10)-(12) 
since the number of parameters remains the same and only the final form of the 
matrices will be different. As concerns Eq.(ll) we made the choice that leads to 
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the simplest form for the matrix elements of Wn as polynomial functions of sines 
and cosines which enter the parametrization of orthogonal matrices. For example 
instead of Wn = On dn_i Wn-i we could take Wn — On Wn-i, where W„_i is at its 
turn given by a formula like Eq.(ll) and so on. 
Examples. An element A4 E U (4) factors as 

A4 = ci4 O4 d\ Ol 4 Ol d\ 

where = {e''^\e''^\ e''^\ e''^^) 4 = (1, e''^^,e''^^, e''^') dl = {l,l, e''^\ e''^^), df = 
(1, 1, 1, e''^'°) and O4, Ol and O^ are the following matrices 



O4 



I cos 9i 

sin 61 cos 62 
sin 61 sin 62 cos C3 



—sinOi 

cos 61 COS 02 

COS 9i sin 02 cos 6^ 





— sin 02 
cos 02 cos 63 







—sin C3 



V 



sin Ui sin 02 COS cos Ui sin 02 cos cos 02 cos —sina^ 
sin 01 sin 02 sin 0^ cos 0i sin 02 sin 6*3 cos 02 sin 0^ cos 0^ J 



01 = 



( 1 

COS 04 

sin 04 cos 0f, 

\ sin 04 sin 0^ 





— sin 04 
cos 04 cos 05 
cos 04 sin 05 



01 = 



The formula (10) for O4 is 



^10 

1 

cos0(i 

\ sin0Q 






—sin0Q 
cos 0Q J 



\ 



—sin 05 

cos 05 J 

\ 



(\ 

1 









cos0^ 
y sin0^ 



O4 — J3,4 • J2,3 • J\,2 — 
\ / 1 \ 

COS 02 —sin 02 

-sin03 sin 02 cos 02 

cos 03 J \ o 1 y 



/ COS 01 —sin 01 \ 



sin 01 cos 01 



\ 















1 

ly 



The Weyl form of a 2 x 2 unitary matrix is 

W2 = w* d2 W2 = dl* Ol d2 O2 dl = 

e*"^! cos^ + e''^2 sin^ cos sin e''^^ (e''^^ 



cos sin 9 e-''^^ (e^'^i - e''^^ ) e''^^ cos^ 9 + e''^^ sin^0 
where = (e^^S e^^^)^ ^1 = giv's) ^-^^ = jj^ ^-^^ u jg the matrix (1). 
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4 Concluding remarks 



In this paper we proposed a new factorization of unitary matrices which can be 
useful in many domains of mathematical and theoretical physics. For example 
the construction of coherent states for SU{n) can be developed similarly to that 
for the SU{3) group given in [0]. The SU{n) matrices are obtained from the 
above parametrization by imposing the condition ipi = upon the 

phases and let us denote by the same letter An an arbitrary element of SU (n) . 
We consider the first n — 1 columns of An as vectors Vk whose entries are given 
by 

{vk)i = Ai^k , i = 1, . . . ,n, k = 1, . . . ,n - 1 

and consider a set of n{n — 1) annihilation operators that we view as the compo- 
nents of (n — 1) annihilation vector operators 

(afc)i = Oifc, 2 = 1, . . . , n, A; = 1, . . . , n - 1 

with the usual commutation relations 

[dik , o.ji] = 0, [ojfc , Oji] = 6ij 6ki 

The key element in defining coherent states is the generating function 

n-1 



■ ■ ■,Vn^i >ni,...,n„_i= a/^i ! • • • n„_i! cxp ^(a]',fi) |0 > 



where ni, . . . , n„_i are the integers that index the representation, |0 > is the vac- 
uum vector and we used the standard notation f for the adjoint of the annihilation 
operator. 

Another problem could be the finding of Laplace-Beltrami operators on uni- 
tary groups that is an old problem |jT^. The Laplace-Beltrami operator on 6*2^-1 
can be written easily as 



A = V - —(cos e,. sin^^''-^^-^e^—) + 



1 ^ 



1 



k=l 



sin^Oi . . . sirPOk-i cos'^Ok dip\ sin'^Oi . . . sin^On-i dip^ 



where we used the parametrisation f„ = (e'^icos6'i, . . . , e'^^sinOi . . . sin6n-i) G 
'S'2n-i With factorization (10) the mathematical tractability problem for SU{'i) 
|l^ and other unitary groups can be resolved. A complete treatment of such 



problems will be given elsewhere. 
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Our proposal is largely based on several simple ideas suggested earlier by many 
people and we have written the factorization in the simplest possible way. We 
suggest it to become the standard one since any other existing parametrization 
can be brought to this form by multiplying with permutation matrices like Pij 
and/or diagonal matrices with entries containing phases. 
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